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ABSTRACT 

Aims. In order to study the acceleration and propagation of bremsstrahlung-producing electrons in solar flares, we 
analyze the evolution of the flare loop size with respect to energy at a variety of times. A GOES M3.7 loop-structured 
flare starting around 23:55 on 2002 April 14 is studied in detail using Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI) observations. 

Methods. We construct photon and mean-electron-flux maps in 2-keV energy bins by processing observationally-deduced 
photon and electron visibilities, respectively, through several image-processing methods: a visibility-based forward-fit 
(FWD) algorithm, a maximum entropy (MEM) procedure and the uv-smooth (UVS) approach. We estimate the sizes 
of elongated flares (i.e., the length and width of flaring loops) by calculating the second normalized moments of the 
intensity in any given map. Employing a collisional model with an extended acceleration region, we fit the loop lengths 
as a function of energy in both the photon and electron domains. 

Results. The resulting fitting parameters allow us to estimate the extent of the acceleration region which is between 
~ 13 arcsec and ~ 19 arcsec. Both forward-fit and uv-smooth algorithms provide substantially similar results with a 
systematically better fit in the electron domain. 

Conclusions. The consistency of the estimates from these methods provides strong support that the model can reliably 
determine geometric parameters of the acceleration region. The acceleration region is estimated to be a substantial 
fraction (~ 1/2) of the loop extent, indicating that this dense flaring loop incorporates both acceleration and transport 
of electrons, with concurrent thick-target bremsstrahlung emission. 

Key words. Sun: flares — Acceleration of particles 



1. Introduction 

Solar flares are known to produce large quantities of accel- 
erated particles, in particular electrons in the deka-keV to 
deci-MeV range. However, the location and physical proper- 
ties of the acceleration region are yet to be well constrained. 
An intrinsic complication is that the radiation produced by 
energetic particles emanates not only from the accelera- 
tion region itself, but also from other locations in the flare 
into which the accelerated parti cles propagat e. Indeed, the 
oft-used "thick-target" model (Br own! [l97lft exploits this 
very complication by deriving properties of the hard X-ray 
emission that are completely independent of the location, 
extent, or physical properties of the acceleration region. 
Hence, determination of the properties of the acceleration 
region from spatially-integrated observations of flare emis- 
sion is not straightforwardly possible. The reader is referred 
to rece nt reviews on electro n properties inferred from hard 
X-rays ([Kontar et al. 2011a ) and the ir implications for elec- 
tron transport ([Holman et al. 2011). 

With the availability of high-quality hard X-ray imag- 
i ng spectroscop y data from the RHESSI instrument 
(iLin et al.l l2002t), the s ituation has much improved (see, 
e.g., lEmslie et all 12003ft . Higher energy electrons are able 
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to propagate further from the acceleration region and 
hence produce hard X-ray emission over a greater sp atial 
extent than i n lower -energy bands. IXu et al] ([2008ft and 
iKontar et al.l (|2011bft analyzed a set of events character- 
ized by simple coronal flare loop sources located near the 
solar limb. In order to determine the spatial properties of 
the flare loops they fitted the RHESSI visibilities with the 
geometric parameters of the loops and determined the size 
of the acceleration regions by fitting the source extents as a 
function of the photon energy with a collisional acceleration 
and propagation model. 

The present paper extends this kind of analysis. 
Specifically, for the simple coronal loop event observed by 
RHESSI on 2002 April 14, we study the variation of source 
extent with energy not only in the photon energy domain, 
but also, for the first time, in the electron domain, using the 
procedure for gener ating weighted mea n electron flux maps 
first enunciated by iPiana et al.l ([2007ft . This extension of 
the model to the electron domain not only admits a sim- 
pler description of the source size with energy E but also 
allows us to exploit the "rectangular" nature of the spectral 
inversion process (in bremsstrahlung, photons with a given 
energy e are produced by electrons with all higher energies 
E > e), making it possible to reconstruct electron maps at 
energies higher than the maximum photon energy observed 
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(see lKontar et al.l l2QQ4). Moreover, the regularization algo- 
rithm used to invert count visibilities into electron visibil- 
ities introduces a natural smoothing over energy E, which 
leads to a smoother behavior of source size with energy and 
so a more reliable estimate of pertinent parameters such as 
the acceleration region length. 

For a given time interval, we use both photon and elec- 
tron maps to examine the form of the variation of loop 
length with energy. The analysis employs three different 
visibility-based imaging algorithms for both photon an d 
electron maps: visibi lity forward-fit ([S chmahl et al. 2 0071), 
maximum entropy (iCornwell & Evansl Il985t | Bong et al.l 
12006( 1. and uv-smooth ([Massone et al.l 12009( 1. In the 
visibility- forward- fit procedure, the loop sizes are deter- 
mined as model parameters. For the other methods, we use 
the standard deviation - square root of the second normal- 
ized moment - of the intensity map as a measure of the 
loop half-lengt h. We then use th e "tenuous acceleration re- 
gion model of IXu et al.l ((2008) to derive the longitudinal 
and lateral extents of the acceleration site. In Section [2j 
we describe the inversion algorithms used to derive elec- 
tron flux visibilities and the imaging techniques employed 
to create the corresponding electron maps for a given time 
interval and electron energy range. In Section [3j we fit the 
variation of loop size with electr on energy E to a simple 
parametric model (|Xu et al.ll2008f ) in order to determine the 
longitudinal and lateral extents of the acceleration region. 
In Section |4j we compare the values obtained through dif- 
ferent imaging techniques and from different map domains 
(photon and electron). 

The inferred length of the acceleration region (~ 15 arc- 
sec) is approximately half the total length of the loop. 
This suggests that the standard model of solar flares where 
electrons are initially accelera ted at a reconnection s ite 
near/above the loop top (e.g., iKopp fc Pneuma nl 119761) is 
not appropriate for certain types of flares. In such flares, the 
acceleration instead takes place over a large region inside 
the flare loop. 

2. RHESSI Visibilities and Imaging Processes 

Solar flare hard X-ray emission is principally produced 
by ac celerated elec trons through the bremsstrahlung pro- 
cess (|Brownl (l97lf ). The relation between the mean elec- 
tron flux spectrum F(x,y;E) (electrons cm -2 s -1 keV -1 ), 
averaged over the line-of-sight direction z through the 
point (x, y) in a target source, and the corresponding 
bremsstrahlung hard X-ray spectrum I(x, y\ e) (photons 
cm -2 s -1 keV -1 arcsec - 2 ) emitted per un it area of the 
source, can be written as (|Piana et al.ll2007l ) 

a 2 r°° — 

V\ e) = J e _ N(x, y) F(x, y; E) Q(e, E) dE, (1) 

where a = 7.25 x 10 7 cm arcsec -1 , representing the extent 
of a one arcsec source at a distance R = 1 AU; Q(e,E) 
(cm 2 keV -1 ) is the bremsstrahlung cross-section for emis- 
sion of a photon at energy e; and N(x,y) (cm -2 ) is the 
column density along the line-of-sight direction. 

By recording the temporal modulation of the de- 
tected flux passing t hrough se t s of rotating absorb- 
ing grids, RHESSI (|Lin et al.l l2002f ) encodes imag- 
ing information in terms of a set of spatial Fourier 



components of the source, termed visibilities, dis- 
tributed over nine circles in the spatial frequency (u,v) 
plane. We define the count visibility spectrum V(u,v;q) 
(counts cm -2 s -1 keV -1 ) as the two-dimensional spa- 
tial Fourier transform of the count spectrum image 
J(x,y;q) (counts cm -2 s -1 keV -1 arcsec -2 ). Similarly, 
the electron flux visibility spectrum W(u,v;E) (elec- 
trons cm -2 s -1 keV -1 ) represents the two-dimensional spa- 
tial Fourier transform of the line-of-sight-column-density- 
weighted mean electron flux image N(x, y) F(x, y; E) (elec- 
trons cm -4 s -1 keV -1 ). 

The relation between the observed count visibility spec- 
trum V(u, v; q) and the e lectron flux visibility spectrum 
W(u, v\ E) is ([Piana et al1l2007l ): 

1 f°° 

V ( u > v >Q)=4^]p J W(u,v;E)K( qj E)dE. (2) 
Here the kernel K(q, E) satisfies 

K(q,E)dq= / D(q,e)Q(e,E)de, (3) 

J e=q 

where D(q, e) is the detector response matrix. 

In order to invert Equation (j2j) to obtain the electron 
flux visibility spectrum W(u, v; E) from the observed count 
visibility spectrum V(u, v; g), we employed a Tikhonov reg- 
ularization technique, which has been proven to be a ro- 
bust and effective inversion method that results in visibil- 
ities (and so images) that vary smoothly with electron en- 
ergy E. Then, from either photon or electron visibility sets, 
we produce the corresponding images using the visibility- 
forward- fit, maximum entropy and the uv-smooth interpo- 
lation/extrapolation methods, as described in Sections I2TT1 
through 1231 respectively. 

2.1. Visibility-based Forward-fit Algorithm 

The visibility-based f orward-fit (FWD) imaging algorithm 
([Schmahl et al.l 120071 ) assumes a parametric source form 
and determines the values of the model parameters that 
result in the best fit to the visibility data. This method 
provides not only quantitative values of the parameters 
but also their uncertainties. Although the applicability of 
the FWD algorithm rapidly deteriorates for complex flare 
morphologies, because of the relatively large number of pa- 
rameters required to characterize the source structure ad- 
equately, the FWD approach is rather effective for sources 
with a relatively simple structure such as the one we study 
here. 

The FWD routines embedded in the RHESSI 
Solar SoftWare (SSW) provide four simple parametric 
source geometries: a circular-Gaussian-distributed sin- 
gle source, multiple Gaussian sources, an elliptical- 
Gaussian source, and a curved-el l iptica l- Gaussian loop 
(jHurford et all 120021 : ISchmahl et aDl2QQ7h . Since the flare 
we are analyzing has a simple loop-structured geometry, we 
adopt the curved elliptical Gaussian form. One of the pa- 
rameters determined by this routine is the FWHM, which 
for a Gaussian profile, is related to the standard devi- 
ation a by exp [-(FWHM/2) 2 /2a 2 ] = 1/2, i.e., a = 
FWHM/V8\K2. 
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We applied this method to visibility data in both the 
photon and electron domain to obtain the standard devi- 
ations cr(e) and <j(E) in the photon and electron domains, 
respectively. Since photons of energy e are produced by 
electron of all energies E > e, these two standard devi- 
ations are not the same; indeed, one expects in general 
a(e = E) > a(E) (see equations [TO] and [T2]V 

2.2. Maximum Entropy Method 

The basis of the Maximum Entropy Method is to maxi- 
mize the information entropy while minimizing the x 2 of 
fit and maintaining the c orrect value of th e total flux. In 
the MEM-NJIT method (|Bong et al.l l2QQ6h . the visibility 
amplitudes V are used to calculate the overall flux of the 
map. The method implements a statistical regularization 
method where the functional 



spectively, can be found by considering the pertinent second 
normalized moments: 



H 



J = H-a X 2 



1=1 



(4) 
(5) 

(6) 



is minimized by means of an iterative scheme. 

In these equations, / is the X-ray image made of n P i X 
pixels, V and V are the observed and predicted visibilities, 
respectively, and the a denote the standard deviations asso- 
ciated with each visibility. The regularization parameter a 
is obtained by means of optimiza tion techniques. T he MEM 
algorithm implemented in SSW ([Bong et al.ll2006l ) provides 
reliable reconstructions in the case of compact events, al- 
though the method often super-resolves the sources and can 
present convergence problems. 



2.3. Fourier-based UV-Smooth Imaging 

An alternative visibility-based imaging method, termed the 
uv -smooth (UVS) algori thm, has been recently developed 
bv lMassone et all ([20091 ) and is available on the SSW tree. 
This method first interpolates the sparsely distributed vis- 
ibilities to generate a smooth continuum of Fourier com- 
ponents in the spatial-frequency v) plane. Then it per- 
forms an FFT-based constrained iterative algorithm to ob- 
tain out-of-band extrapolations. The method has proven 
to reproduce realistic forms of the source structure with 
a high degree of accuracy, fidelity, robustness, and com- 
putational efficiency. Although the method may introduce 
artifacts when applied to source configurations character- 
ized by distant footpoints, uv-smooth is very accurate when 
reconstructing relatively localized extended sources such as 
those considered here. 

Unlike images produced by the FWD algorithm, the 
UVS and MEM procedures do not straightforwardly pro- 
vide quantitative information on the uncertaintiesjn the de- 
termined source extents. Given the electron flux F(x,y;E) 
in a two-dimensional flux image at any given energy E, the 
location of the source can be estimated by calculating the 
first normalized moment of the intensity. The length and 
width of a source, in the photon and electron domains re- 



a(9;E) 



f™sH(s,0;e)ds 
"V S™I{s,6;e)ds 

' s 2 N(s)F(s,0;E) ds 
/ °° N(s)F(s,6;E)ds ' 



(7) 

(8) 



where (s, 6) are polar coordinates in the plane of the sky, 
relative to an origin that we define as the location of the 
maximum flux intensity. The integrals in Equation ([7]) are 
computed numerically for a variety of 6 values. The max- 
imum (minimum) values of a (9) can be identified as the 
length (width) of the source. In order to provide quantita- 
tive uncertainties on these values, we applied a Monte Carlo 
approach in which random noise is added to the visibilities 
and the resulting images recomputed and reanalyzed. Ten 
realizations of the source visibilities were used; the stan- 
dard deviations of these ten results were taken to be the la 
errors of the loop length and width. 

3. Application to A Loop-Structured Flare 



In the standard flare model, thick-target "footpoints" 
are considered to represent the dominant locations of hard 
X-ray emission because the coronal magnetic loops through 
which the energetic particles propagate are generally not 
dense enough to stop electrons via Coulomb collisions. 
Due to the high density of the chromosphere, the hard X- 
ray structure of footpoint sources typically extends over a 
very small spatial extent, so that observations with avail- 
able spatial resolutions cannot directly determine details of 
the particle acceleration and propagation processes in the 
bremsstrahlung-emitting region. 

However, RHESSI has revealed a new class of flares in 
which the hard X-r ay emission is predominantly from the 
coronal loop itself (iKrucker et al1l2008l ; IVeronig fe Brownl 
l2004t ISui et al"1l2004l ). For such sources, the corona is not 
only the location of particle acceleration, but also dense 
enough to act as a thick target, stopping the accelerated 
electrons before they can penetrate to the chromosphere. 
Close investigations of these flares can provide direct infor- 
mation on the electron acceleration and propagation pro- 
cesses in the bremsstrahlung-emitting region. 

One of the most closely studied events from this class 
is the "midnight" flare of 2002 April 14 (ISui et all 120041: 
Veronig fc Brownl 120041 : iBone et al.l 120071 : IXu etafl 120081 : 
Kontar et al.ll2011bl ). Figure [Qi shows the RHESSI count 
rate profiles for this event in five different energy channels. 
For its entire duration this flare consists of a simple loop 
structure viewed "side-on" near the solar limb, as shown in 
Figure [T^. The c oronal loop has a densit y sufficiently high 
(n ~ 10 11 cm -3 ; IVeronig fc Brownl |2004( ) to stop electrons 
up to ~ 30 keV over the observed loop length of 10 9 cm. It 
has also been suggested that electrons must be continuously 
accelerated during this event because for such a dense loop 
20 keV particles hav e collisional lifetimes less than 0.1 s 
([Kontar et al.ll2011bl ). 
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Fig. 1: [a]: Time profile in different RHESSI energy ranges 
of the flare starting around 23:55 UT on 2002 April 14. 
The vertical lines indicate the five time intervals for which 
the spectral fitting and visibility-based imaging were per- 
formed, [b]: The spatially integrated RHESSI X-ray spec- 
trum fitted with an isothermal plus thick-target nonthermal 
model for the third time interval (00:06 - 00:08 UT). [c]: An 
example of a photon flux map (for the third time interval 
and energy bin 18-2 keV) obtained by the Clean method 
(jHurford et aLll2QQ2h - 
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Fig. 2: Mean electron flux maps for the time interval 00:06 
- 00:08 UT, obtained through visibility- forward- fit (left ), 
Maximum Entropy (MEM-NJIT) (middle) and uv-smooth 
(right) procedures applied to electron visibilities. Energy 
bins are from 12-14 keV (top) to 28-30 keV (bottom). The 
range for x-axis is [744, 824] arcsec and the range for y-axis 
is [342, 422] arcsec. 



We have analyzed both the photon flux maps and 
the mean electron flux maps for five different time inter- 
vals (00:02-00:04, 00:04-00:06, 00:06-00:08, 00:08-00:10, and 
00:10-00:12 UT on 2002 April 15) and ten 2-keV energy 
bins from 10-12 keV up to 28-30 keV. Visibilities sampled 
by the front detector s of collimators 3 to 9 are used for th e 
map reconstruction (|Xu et al.l 120081 : iMassone et al.l l2009f ) . 
Figure [2] shows the mean electron flux maps at time inter- 
val 00:06 - 00:08 UT, obtained through FWD, MEM-NJIT 
and UVS procedures, all using the electron visibilities ob- 
tained thro ugh spectral invers ion of the observed count visi- 
bility data ([Piana et al.ll2007l ). All three methods reveal the 
flare structure as a single loo p. It should be noted that in 
general as has been shown by IMassone et all ([200 9), using 
simulated data consistent with plausible astrophysical con- 
ditions, MEM-NJIT tends to underestimate source sizes. 
We therefore use only FWD and UVS algorithms for the 
analysis to follow. Figure [3] shows the UVS mean electron 
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where L(e) (arcsec) is the loop length at photon energy e 
(keV); I/o (arcsec) is the extent of the acceleration region, 
and a (arcsec keV -2 ) is a parameter inversely proportional 
to the plasma density n in the flare loop: 



1 



a 



(5-2) 



Kn (<5-3)(<5-4)' 



(10) 



Here K = 27re 4 A (e being the electronic charge and A being 
the Coulomb logarith m) and 5 is the spectral index of the 
injected electron flux (|Xu et al.ll2008l ). This model assumes 
that electrons are accelerated within a tenuous region ex- 
tending from —Lo/2 to Lo/2 and are injected into an exter- 
nal region where the loop density is sufficiently high that 
propogating electrons may lose energy through Coulomb 
collisions and produce hard X-ray emissions Q. Moreover, 
electrons with higher energies can propagate further and 
hence produce photon emission over a larger spatial ex- 
tent. Employing a similar analysis, we can extend the above 
photon-based model to the electron domain. This yields the 
result: 



L(E) = L + PE 2 , 



(11) 



where L(E) (arcsec) is the electron loop length at electron 
energy E (keV) and 



00:04-00:06, 



00:08-00:10, 28-30keV 






(3 = 



Kn (6 - 3) ' 



(12) 



Since the emission at photon energy e is a weighted sum 
of electron flux at energies E > e, the overall loop extents 
in the photon domain (Figure H]) are generally largei0 than 
those in the electron domain (Figure [SJ. 



Fig. 3: Mean electron flux maps obtained through uv- 
smooth algorithm (Section 12. 3p for time intervals 00:04 
- 00:06 (left), 00:08 - 00:10 (middle), and 00:10 - 00:12 
(right) UT. Energy bins are from 12-14 keV (top) to 28- 
30 keV (bottom). The ranges for x-axis and y-axis are the 
same as in Fig. [2j 

flux maps in five different energy bins (every second 2-keV 
energy bin from 10 to 30 keV) and three different time in- 
tervals. It can be seen that the loop length continuously 
grows with energy. However, with time, the size differences 
between low and high energies become smaller, indicating 
(see below) a loop density that increases with time. 

For the FWD imaging process we adopt a curved- 
elliptical-Gaussian model, for which the loo p geometry is 
modeled by seven parameters (|Xu et al.ll2008f ). One of these 
seven fitting parameters, i.e., the source extent parallel to 
the curved arc and termed the length of the loop L, is par- 
ticularly important to this study. For the UVS maps, the 
loop length is estimated from calculating the second nor- 
malized moment of the flux intensity using Equation (0. 
Figures |4] and [5] show the lengths of the source as a func- 
tion of photon energy e and electron energy E for five given 
time intervals throughout the flare impulsive phase. The 
loop lengths clearly grow wi th energy. In the "tenuous ac- 
celeration region" model of IXu et al.l ([20081 ), the form of 
L(e) is 

L(e)=L + ae 2 , (9) 



Figure [TJd shows the spatially integrated RHESSI X- 
ray spectrum at the third time interval (00:06 - 00:08 
U T) obtained using all front detectors excluding 2 and 
7 ([Smith et al.l l2002f ). By fitting the spectrum with an 
isothermal pl us collisional thick-target nonthermal model 
([Brownl Il97lf ) , we have found that the transition energy 
between the themal and nonthermal components is about 
15 keV during this event. Therefore, we fit the acceleration 
model in both the photon and electron domains (equations 
[9] and [11]) starting from E = 14 keV. Spectral inversion 
of the photon-based data allows reconstruction of electron 
maps at energies (up to 30 keV) higher than the photon 
energies (up to 26 keV) observed. 

The loop lengths and fits based on both FWD and UVS 
methods are presented in Figure [4] (photon domain) and 
Figure [5] (electron domain). The best-fit acceleration re- 
gion lengths Lo, and their statistical uncertainties (given 

1 Note that a more correct "dense acceleration region" form 
for L(e), which in corporates emis sion from the acceleration re- 
gion itself, exists (|Xu et al.|[2Q Q8f) . The pertinent form of L(E) 
is more difficult to use in a best-fit analysis; however, this 
model yields results for Lq and n that are very similar to the 
more straightforward-to-apply "tenuous acceleration region" re- 
sult ©. 

2 From Equations [9] and [11] , the ratio of the propagation 
lengths is ae 2 //3E 2 = ([8 - 2}/[S - 4)]) x (e/E) 2 . 
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Fig. 4: Flare loop lengths for each 2-keV photon energy bin, 
for five different time intervals throughout the flare impul- 
sive phase, [a]: results from forward- fit algorithm; [b]: re- 
sults from uv-smooth method. The solid curves represent 
the fits of the model described in Equation (j9]). The y 2 °f 
the fittings are also shown. 
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Fig. 5: Flare loop lengths for each 2-keV electron energy 
bin, for the same five time intervals as in Figure HI [a]: 
results from forward- fit algorithm; [b]: results from uv- 
smooth method. The solid curves represent the fits of the 
model described in Equation ([TT]) . The x 2 of the fittings 
are also shown. 



by Monte Carlo simulations), are shown in Table [TJ For 
each time interval, we have also estimated the loop width 
W averaged over different energies. Assuming W is the ap- 
proximate extent of the acceleration region across the mag- 
netic loop and the loop is essentially a cylindrical column, 
we obtain the volume of the acceleration region to be 

V = ttW 2 L /4. (13) 



The values of W and Vb, and their statistical uncertain- 
ties, are also shown in Table HJ From these results, it can 
be deduced that: 

- the models describe the data accurately in all contexts: 
both the photon and electron source lengths are well-fit 
by a quadratic form, with similar values of Lq in both 
FWD and UVS results; 



— the accelerat ion region length Lq is comparable to that 
obtained by IXu et al] ([200 8) (~ 12.5 arcsec) using a 
forward-fit algorithm in the photon domain; the slight 
discrepancy could be due to the choice of different algo- 
rithms as well as time and energy ranges; 

— the extension of the model to the electron domain is 
novel, allowing smoother flux variations with energy due 
to the regularizing constraints involved in the inversion 
process; this in turn validates the photon-based results 
and represents one step closer to the underlying physics. 
The electron maps generally offer a better fit (i.e., the 
X 2 values are smaller), particularly using the FWD pro- 
cedure. 

— the loop lengths appear to show a systematic change 
throughout the event: the energy- length curves are 
steeper at earlier time intervals and flatter at later time 
intervals, especially in the electron domain shown in 
Figure [U According to the model (Equation [12]) this 
may be due to the increasing density as the flare evolves, 
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Table 1: Acceleration region length Lq (Equations [9] and [TT]). loop width W and volume ^(Equation [13]). 



Time 


(hh:mm) 


00:02-00:04 


00:04-00:06 


00:06-00:08 


00:08-00:10 


00:10-00:12 


Photon Map 




L (FWD) 


(arcsec) 


14.3 ± 1.9 


13.3 ± 1.2 


14.1 ±1.2 


15.0 ±1.0 


16.4 ±0.7 


L (UVS) 


(arcsec) 


13.3 ±1.3 


14.3 ± 1.1 


15.2 ±0.9 


15.1 ±1.1 


16.4 ±0.9 


W (FWD) 


(arcsec) 


7.1 ± 1.1 


8.6 ±1.0 


8.1 ±0.8 


7.1 ±0.6 


6.8 ±0.5 


W (UVS) 


(arcsec) 


7.0 ±1.5 


7.5 ±1.3 


7.7 ±1.5 


6.6 ±1.3 


7.3 ±1.1 


V (FWD) 


(100 arcsec 3 ) 


5.7 ±1.9 


7.7 ±1.9 


7.2 ±1.6 


5.9 ±1.1 


6.0 ±0.9 


(UVS) 


(100 arcsec 3 ) 


5.1 ±2.2 


6.3 ±2.2 


7.1 ±2.8 


5.2 ±2.1 


6.9 ±2.1 


Electron Map 




L (FWD) 


(arcsec) 


15.9 ±0.9 


13.3 ±0.8 


16.3 ±0.9 


17.4 ±0.8 


18.3 ±0.6 


L (UVS) 


(arcsec) 


15.7 ±0.7 


14.7 ±0.8 


16.5 ±0.6 


17.0 ±0.9 


18.1 ±0.6 


W (FWD) 


(arcsec) 


5.8 ±1.5 


7.8 ±1.5 


7.4 ± 1.3 


6.5 ±0.8 


6.2 ±0.6 


W (UVS) 


(arcsec) 


7.7 ±1.8 


7.1 ±1.5 


7.4 ±1.7 


7.0 ±1.6 


7.6 ±1.6 


V (FWD) 


(100 arcsec 3 ) 


4.2 ±2.2 


6.4 ±2.4 


6.9 ±2.6 


5.7±1.5 


5.5 ±1.0 


Vb (UVS) 


(100 arcsec 3 ) 


7.3 ±3.4 


5.8 ±2.5 


7.1 ±3.2 


6.5 ±3.1 


8.1 ±3.4 



given that the electron spectral index 5 does not show 
a systematic change with time 0. 



4. Interpretation of the Results 

It is often postulated that electrons are accelerated at a re- 
connection current sheet and that the acceleration process 
is decoupled from the electron transport that follows the 
reconnection (|Aschwandenlll998l ). Indeed, many solar flares 
dem onstrate a clear loop-top-pl us-footpoints structure (see, 
e.g.. iBattaglia fc Kontarll201lf ). Howeve r, this flare has a 
loop density as high as n ^ 10 11 cm -3 ([Veronig fc Brownl 
\2004 ). Therefore, most of the energetic electrons would 
be stopped through Coulomb collisions near the looptop 
([Kontar et al.l l2011bf ) , and so the emission region would 
be confined to a limited extent less than the loop length 
(~ 30"). The spatially resolved observations presented here 
suggest the presence of an extended acceleration site inside 
the loop. Further, given that the acceleration region (<~ 15") 
is estimated to be a substantial fraction of the total length 
of the flare loop (~ 30"; see Figures H] and [5]) and that the 
hard X-ray emission is produced in the entire loop instead 
of only at the footpoints, it is clear that the dense flare loop 
incorporates both acceleration and transport of electrons, 
with concurrent thick-target bremsstrahlung emission over 
the entire length of the source. 

Such a scenario is consistent with the presence of 
an enhanced leve l of MHD fluctuations inside the loop 
(|Bian et al.l l201ll ) . It is possible that for such events re- 
connection and loop-top injection may not correspond to 
the dominant energy dissipation and particle acceleration 
process, which instead may proceed inside the flare loop it- 
self t hrough stochastic mechanisms (e.g., iPetrosian fc Liul 
l2004f ). Our result shows that a considerable proportion of 
the flare loop is involved in the pr ocess of acceleratio n, thus 
avoiding the "number problem" ([Brown et al.1120091 ). 



3 The spectral index of the injected electron flux S is obtained 
from RHESSI X-ray spectra fitting using an isothermal plus 
thick-target nonthermal model. Through the five studied time 
intervals, 5 has values as [7.6, 7.3, 7.9, 8.6, 7.8] 



5. Conclusions 

We have analyzed an extended coronal hard X-ray source 
in the event of 2002 April 14 with photon maps con- 
structed from count visibilities (two-dimensional spatial 
Fourier transforms of the source geometry). Furthermore, 
using a regularized spectral inversion technique generat- 
ing electron visibilities, we have studied the dependence of 
the electron flux image on electron energy E. The source 
lengths derived from both photon and electron maps gen- 
erally grow quadratically with energy, in agreement with 
a col lisional model i nvolving an extended acceleration re - 
gion (|Xu et al.ll2003lKontar et al.ll2011bHBian et al.ll201lh . 
Fitting this model allows estimation of the length and vol- 
ume of the acceleration region: ~ 15 arcsec and ~ 600 
arcsec 3 , respectively. We compare the results obtained by 
different algorithms (FWD and UVS) and in both photon 
and electron visibility domains. The plausible and consis- 
tent estimates of the acceleration-region lengths and widths 
strongly suggest that the proposed model is reliable. The 
systematic change of the behavior of loop length with time 
indicates that the loop density may increase th roughout the 
event , due to chromospheric "evaporation" ([Acton et al.l 
1982). The inclusion of a thermal component to the cur- 
rent model is required for a comprehensive description of 
the chromospheric evaporation process. 

Detailed and statistical studies of more cases of dense- 
loop flares will be carried out in future work, in which we 
will also investigate other properties of the acceleration re- 
gion, e.g., the loop density and its evolution with time. 
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